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Abstract 

A general method for computing 0(a 2 ) and higher-order next-to- 
leading photonic corrections is presented and applied to the precision 
calculation of the small-angle Bhabha scattering cross section in the 
phase-space region of interest for the luminosity measurement at LEP. 
The formulation is based on a proper matching of exact 0(a) results 
with higher-order corrections in the Structure Function formalism. The 
results of the approach are analytically compared with theoretical cal- 
culations, both for s- and t-channel processes, available for simple Event 
Selections. Numerical predictions for realistic Event Selections are also 
provided and critically compared with the ones existing in the literature. 
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1. Introduction 



Small-angle Bhabha (SABH) Scattering is used at e + e~ colliders to measure 
the accelerator luminosity. Since the experimental error at present reached by 
the LEP Collaborations is better than 0.1% |]T|, in order to exploit such an 
experimental achievement it is mandatory to provide theoretical predictions 
for the cross section in the luminosity region at the same level of accuracy, or 
even better. 

During 1995, within the Workshop on "Physics at LEP2" at CERN, the 
Working Group "Event Generators for Bhabha Scattering" took place ||. The 
major task of the Working Group was to make an inventory of all the avail- 
able Monte Carlo (MC) and non-MC codes for the computation of the SABH 
cross section, with the aim of reaching a deeper understanding of the SABH 
process. The main result of the Working Group was a substantial reduction 
of the theoretical error on the SABH cross section |2j, ||| from 0.16% to 0.11%, 
but further improvements are demanded. In particular, at the present stage 
the dominant part of the theoretical error comes from the 0(a 2 L) photonic 
corrections, where L = ln(Q 2 /m 2 ) is the usual collinear logarithm, which for 
the time being are not fully under control for a realistic experimental Event 
Selection (ES). 

The main tool by means of which such a result has been achieved, was 
the critical comparison of all the available theoretical formulations and the 
corresponding numerical codes [Q, ^J, ^, [7|, and the understanding of their 
differences, where present. Given the relevance of this subject, we considered 
important to improve our theoretical formulation presented in 0, with the 
goal of providing a contribution to a further reduction of the theoretical error. 
The aim of the present note is to describe the recent theoretical developments 
of our approach to the understanding of the SABH process, with particular em- 
phasis on 0(a 2 L) and higher-order next-to- leading photonic corrections. The 
results of the approach are analytically compared with theoretical calculations, 
both for s- and t-channel processes, available for simple ES's. Numerical pre- 
dictions for realistic ES's are also provided and critically compared with the 
ones existing in the literature. 

From now on, we will focus our attention on QED radiative corrections 
to the dominant part of the cross section in the small-angle regime, namely 
the cross section for t-channel photon exchange. The corrections due to the 
other non-dominant contributions, namely s-channel diagrams, j-Z interfer- 
ences and so on, can be accounted for at the leading logarithmic (LL) level, 
as shown in 0. We will start considering Bare ES's. They are absolutely 
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unrealistic from the point of view of luminometry at LEP, so their interest is, 
in some sense, academical. In spite of this, they however represent a useful 
benching situation, since the only available calculation complete at the level 
of 0(a 2 L) corrections is conceived for such ES's. The results concerning 
Calorimetric ES's, the truly realistic ones, will be derived as special cases of 
the ones obtained for Bare ES's: as a first step, final-state radiation will be 
switched off in the LL part of the result, consistently with the fact that, for 
an inclusive measurement, final-state mass "singularities" are absent; next, it 
will be taken into account by means of a proper definition of the final-state 
radiation factor. 

2. Theoretical formulation 

In order to settle down the basic notations, let us recall the form of the 
corrected cross section in the LL approximation within the Structure Function 
(SF) formalism |§: 

a LL = J dx 1 dx 2 dy 1 dy 2 J dI cm ,Dp.(xi)Dp.(x 2 ) 

(I cm (x 1 ,x 2 ,y 1 ,y 2 ))D f3j (y 1 )Df 3f (y 2 ). (1) 



dl 



cm 



In eq. (P, Dp^xip) are the usual electron and positron structure functions, 
describing initial-state radiation, as can be found in |7]]; (2/1,2) are the anal- 
ogous ones for final-state radiation. In general, /3, and (3f are the radiation 
factors for initial- and final-state radiation, respectively; in the case of a Bare 
ES, one has /3j = f3f = 2(a/7r) [ln(— t/m 2 ) — 1]. The case of a Calorimetric ES 
will be discussed later on. da/dl cm is the density to be QED corrected, in the 
centre of mass (cm) reference frame after initial state radiation. I{x\, x 2 , yx, y 2 ) 
represents a convenient set of independent variables, able to describe the con- 
figuration of the event corresponding to the kernel cross section, in presence 
of initial- and final-state radiation. It is understood that the event generated 
in the cm frame after initial state radiation is boosted along the beam-axis 
according to the amount of radiation lost longitudinally by the incoming par- 
ticles, then modified by the effect of final-state radiation, and then accepted 
or rejected according to a given ES in the laboratory frame. 



The 0(a) content of eq. (|T|) can be written as 



a 



(a) 
LL 




dc 

{I cm {x, 1, 1, 1)) 



diem 



3 



da 



dlcrr 



(/ cm (l,x,l,l))-2-^-(/ cm (l, 1,1,1)) 

(J/J. c 1 



da 

dlrrr, 



da 



(1, 1, X, 1)) + — — (I cm (l, 1, 1, x)) 



die 



-2- 



da 



dL 



(7^(1,1,1,1)) 



(2) 



obtained by inserting in eq. (]1|) the distributional expansion of the structure 
functions and dropping the spurious higher-order terms. It is worth noting, 
for future convenience, that the LL cross sections, both the complete and the 
(9(a)-truncated ones, are functionals of the kernel cross section. In particular, 
they define algorithms that, once applied to a generic kernel, provide the whole 
ensemble of the LL corrections to all orders and the 0(a) LL ones, respectively. 



So, when appropriate, we will refer to them as Y,ll da/dl 
respectively. 



nr Y (a) 



da/dl 



The matching between the all-orders leading- log cross section, cr^i, given 
by the convolution of structure functions with kernel (Born) cross sections, 
and the exact 0(a) one is realized according to the following general recipe: 
the 0(a) content of the leading-log cross section is extracted by employing 
the O(a) expansions of the structure functions, thereby yielding a^. By 



denoting with <jpj. the exact 0(a) correction to the cross section |10|, |TTJ, the 
fully corrected cross section can be written as 

<?A = °LL ~ (Tl£ + a PT ■ ( 3 ) 



Equation (|3|), which is in the additive form, is by construction exact at 
O (a) and includes higher order corrections as taken into account by a^. A 
factorized form can also be supplied. It has the same 0(a) content but also 
leads to the so-called classical limit, according to which the cross section must 
vanish in the absence of photonic radiation. It reads 

(Tf = (1 + C^ L )a LL , 



r («) 

Tprp 



a 



(a) 
LL 



_(«) 
a NL 



(4) 

ctq being the Born cross section; C§ L contains the non-log part of the 0(a) 



cross section, represented by a 



(a) 
NL- 



In order to clarify the physical content of eqs. @ and @, it is worth 
considering the perturbative expansions of their ingredients, which can be 
written as follows: 
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„ — rr I n- (o) 4- n- ( ° 2) 4- ^r (a3) 4- 
C>LL - ^0 + O'LL + &LL + a LL +■■■> 

, (<*) 

ctll| q = (To + cr LL ; 

CTpr = O"0 + OpT- (5) 



The dominant part of the difference between the additive cross section and 
the factorized one can be read off the perturbative expansions above: 



*F = *A + Cg L <r% + 0(o?Ce). (6) 

The cross-section a L °£ describes the LL universal part of the 0(pt) corrections, 
and C^l the non-log (process dependent) part of the 0(a) corrections. There- 
fore, the difference between the factorized cross-section op and the additive 
cross-section a a is due to 0(a & £) sub-leading corrections contained in op and 
not in a a- In particular, the term given by C^iPlI gives an approximation of 
the (9(a e £) contributions to the cross section, 

- C H NL atl (7) 

since it is the direct product of the LL 0(aC) cross section times the non- 
logarithmic non-universal 0(a) correction C$ L . A more detailed description 
of the 0(a £ C) correction would require a convolution of the non-universal 
non-logarithmic contribution, given by the difference between the exact O(o) 
cross section minus the LL O(o) one, with real+virtual radiation in the LL 
approximation. This can be achieved by applying the algorithm of eq. (^) in 
the following way: 



a 



(a 2 L) 



<<*) 
J LL 



dai°t) da 



(a)" 
LL 



dVL 



where da/dfl are the differential distributions of the electron scattering angle, 
already integrated over the photonic phase-space. Equations (0) and (§) are, at 
a first sight, very different from one another. On the contrary, eq. @ reduces 
to eq. (0) under the assumption that the ratios between the 0(a) densities 
and the born one are smooth functions of the cm energy and scattering angle, 
which is a priori a reasonable assumption since the 0(a) densities are infrared- 
dominated and hence almost factorized over the born one. Anyway, one of the 
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purposes of the present note is just to check the quality of the approximated 
description of eq. (0) and improve it by means of eq. Some illustrative 
numerical results will be shown later on. 

Equation performs the 0(aC) corrections to the O(a) non-logarithmic 
kernel, and hence provides C(a e £) corrections, taking into account also con- 
volution effects. Along the same way, one can improve it to take into account 
also the 0(a\£\~°°) corrections in the SF approach, by computing, instead 
of (H), its improved version, namely 



a 



dxidx 2 dyidy 2 / dI cm D0 i (x 1 )Dp i (x2)Dp f (y 1 )D/ 3f (y 2 ) 



(J cm (xi, x 2 , 2/1, 2/2)) - -rrK-Wl, !> !> !)) 



(Icm(x 1 ,x 2 ,yi,y 2 )) 



dl c 
da 



(a) 
LL 



dL. 



(J cm (l, 1,1,1)) 



(9) 



Some comments are in order here. If one inserts into eq. the 0{a!) 
contribution of the distributional expansion of the structure function, the result 
is automatically zero, confirming that eq. (||) contains contributions starting 
from 0(a £ C). If one inserts in it the 0(a) contribution of the distributional 
expansion of the structure function, one naturally recovers eq. (U). The higher- 
order contributions in the structure functions produce 0(a K C^~ 00 ) corrections: 
they are complete at the LL level until (9(a 9 £ e ), correct in the soft limit 
for all the higher-order next-to-leading contributions due to the fact that at 
present the structure functions employed take into account only up to C(a e ) 
hard-photon radiation. This, of course, is not a limitation in principle, and 
moreover has no practical influence for the usually adopted experimental cuts. 

The discussion up to now, together with the definition of the various ingre- 
dients, allows the definition of a new, factorized cross section. In fact, it is easy 
to verify that the sum of the additive cross section of eq. (Q) and the higher- 
order next-to-leading correction of eq. (|j) defines the following factorized cross 
section: 



a A + cr y = a p , 



a new = J ( j i x 1 dx 2 dy 1 dy2 J dtt cm D p .(xi)D /3 .(x 2 )Di3 f (y 1 )Dp f (y 2 ) 
-(xi, x 2 , 2/1, 2/2) 1 + C% L (xi, x 2 , 2/1, 2/2), , 



da 

d^lrm 
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(10) 



The analytical identity between a a + cr ( - a ™ L ™ •* on the one hand, and ap ew 
on the other one can be seen as follows. <jp ew is the sum of two integrals. The 
first one, involving only the born-approximation density, is nothing but (Tll- 
The second one, involving the product of the born- approximation density times 
the new C^ L (x\, %2i 2/1, 2/2), coincides with the full convolution of the difference 
between the exact 0(a) density and the LL-approximated one. Let us focus 
the attention on this last integral. By inserting in it the 0(a') distributional 
expansion of the structure functions and performing the integrations, one is 
left with Opj. — which, when summed up to (Jll-, gives exactly a a- The 
residual higher-order contributions to the second integral, by virtue of the fact 
that the structure functions are normalized to unity, recover exactly a ( - a " Ln \ 

Equation fllQ|) reproduces the result of eq. (^) under the assumption that 
C§ L (xi, X2, 2/1, 2/2) is a smooth function of its arguments, so that it can be 
extracted from the convolution integral and the densities appearing in its nu- 
merator and denominator can be replaced by the corresponding integrated 
cross sections. In this sense, eq. ( |I0D is an improvement of eq. (|]), that fully 
takes into account the angle and energy dependence of the non- logarithmic 
component of the exact 0(a) correction, so that eq. (|T0| ) supersedes eq. (f|) 
developed in ref. 0, which the published version of the FORTRAN code SAB- 
SPV U is based on. Moreover, from a technical point of view, the formulation 
described in is worked out under the assumption that a specific cut on 
the final state fermions is applied, namely a cut on the individual energies of 
the fermions. The present formulation allows one, on the contrary, to impose 
arbitrary cuts on the final state fermions. 

The formulation just described is conceived for applications to the SABH 
process. Nonetheless, it is completely general, and can be applied also to 
annihilation processes, once the non-logarithmic 0(a) process-dependent cor- 
rection (Cn L ) is known. 

It is worth mentioning that another source of 0(a e C) corrections is rep- 
resented by the production of additional light pairs; this contribution is at 
present well under control at the level of 0.03% p], [|, and a further improve- 
ment will be necessary only when photonic corrections will be under control 
at the same level. 

So far, only QED corrections to the dominant part of the small-angle 
Bhabha cross section, i.e. the t-channel photon exchange contribution, have 
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been considered. This is not a limitation of the approach: the very same 
algorithm of eq. (|10D can be extended to the full Bhabha cross section in a 
straightforward way. Actually, from the practical point of view, at small scat- 
tering angle (20-60 mrad) it is necessary to take into account the residual 
born-approximation contributions plus their LL corrections, and this can be 
simply obtained by adding the residual born-approximation contributions to 
the integration kernel of eq. fllPf ). Moreover, the only non-QED corrections 
relevant in the luminometry region are given by the photonic vacuum polar- 
ization, which also can be easily taken into account by using the running QED 
coupling constant. 



3. Comparisons and numerical results 

Analytical comparisons - The approach just described insures that the 
0(a G C) photonic contributions coming from an "external" collinear photon 
in association with an "internal" non-collinear one are automatically taken 
into account. This is not in principle the complete set of 0(a € C) corrections: 
in this way, for instance, the truly irreducible two-loop corrections are missed, 
but they can be expected to give small contributions. Therefore, the 0(a G £) 
corrections taken into account by this method represent the bulk of the com- 
plete set. This heuristic argument can be put on a firmer ground by comparing 
the results provided by this formulation with the ones obtained by means of 
complete 0(a e C) calculations, already present in the literature for some spe- 
cific ES's. As a first step, we will compare our results with the ones available 
in analytical form, i.e. in the soft approximation limit; next we will consider 
comparisons beyond the soft approximation, for which only numerical results 
can be compared. 

A first analytical result including 0(a e C) corrections concerns e + e~ anni- 
hilation into pairs, taking into account initial-state QED corrections ]12j . 
In that paper an analytical formula is given, describing QED corrections for 
an ES where the total energy emitted by initial-state photons does not exceed 
AE. If one works out eq. (0), namely by putting in this case f3j = 0, imposing 
the cut condition of |T2j , computing the proper C^ L and truncating the result 



at the (9(a e ), one finds that: 



• the O(a) perturbative result is exactly recovered, by construction; 

• all the infrared-singular terms, namely the ones containing In 2 e and 
hie, where e = AE/E, are exactly recovered at the level of 0(a e £j), 
0(a e Cf) and 0(a e ), where L s = ln(s/m 2 ); 
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the difference between the two results starts at the level of (a/7r) 2 L s 
times a constant; 



in particular, such a difference reads 
5a 



aV 

IX ) 



3C(3)-^C(2) + A 



where 5a is the difference between the cross section of 12 and the cross section 



of eq. (pT0|) . The difference numerically amounts to a relative deviation of 
about 1.7 x 1CT 4 . The residual difference is at 0(a e ) times a constant and is 
numerically irrelevant. 

It is worth noting that results including full 0(a £ £) corrections for the 
small-angle Bhabha scattering cross section are available for an academic trig- 
ger U] or under development for the most general case JT3] . Also in this case it 



is possible to compare analytically the predictions of eq. ( |10D with the results 
shown in M in the soft approximation. The comparison can be worked out 
along the same lines as in the previous case, with the only differences that the 
collinear logarithm is now L t = ln(— t/m 2 ), (3f = = (3 t and the cut condi- 
tion imposed in |J requires that the energy of each photon does not exceed 
AE. The results of the comparison are the same as before up to the 0(a e C u ) 
corrections, namely the difference appears at the level of (a/ix) 2 L t times a 
constant and reads 



5a 



(a 2 Lt) 



1\- 

,7T 



3C(3)-^C(2) + A 



(12) 



where 5a is the difference between the cross section of || in soft approximation 
and the cross section of eq. (p~0|) . This difference numerically amounts to 
a relative deviation of about 2.2 x 10~ 4 , since the overall factor of two is 
compensated by the fact that L t ~ 2/3L s . In this case, also an additional 
difference appears, namely at the level of the infrared- sensitive truly 0(a £ ) 
terms, which reads 

5a 



0-Q 



((a*)) 



a 



71 



4hr e + 8\ne 



(13) 



the residual difference being at 0(a e ) times a constant and hence numerically 
irrelevant. This last difference, being infrared-sensitive, can show up in the 
region of strong cuts. 

Numerical results - Going beyond the soft approximation requires the dis- 
cussion of numerical results. From now on, it is understood that all the nu- 
merical results, apart from the ones obtained by our new formulation, can be 
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found in 0. The formulation described in this paper has been implemented 
in a new, so far unpublished, version of the code SABSPV ||. 

Figure [l] shows the comparison between the results obtained by the present 
formulation of the problem and the results available in the literature |2j for 
a Bare ES. For the details concerning the ES, the reader is referred to ||. 
The FORTRAN codes involved in the comparison are BHAGEN95 [|, BH- 
LUMI 0, NLLBHA @, OLDBIS+LUMLOG [§ and SABSPV 0, g. Again, 
the details concerning the codes and their underlying theoretical formulations 
can be better found in the original literature and in 0. Figure [I] is drawn 
according to the convention adopted in || , namely it shows the relative differ- 
ences between the codes involved and a reference cross section, which in the 
present case is the cross section computed by BHLUMI. As a first comment, 
one can see that our new factorized solution corresponding to eq. (^) differs 
appreciably from the previous factorized one of ref. |7j due to the improved 
treatment of the final state radiation phase space. Both of these solutions, ne- 
glecting convolution effects at the level of 0(a € C) corrections, are at present 
obsolete. By looking at the figure, one can see that all the solutions shown 
group together essentially into two clusters. The first cluster of solutions is 
the one of the "additive" solutions, namely BHAGEN95, OLDBIS+LUMLOG 
and SABSPV in its additive version (see eq. (|3|)). All these solutions have the 
common feature that they miss the 0(a e C) contributions, since they do not 
fill the region of phase space characterized by the emission of one collinear and 
one acollinear photon. Their relative distance from the reference cross section 
is steadily around 0.10-0.15%. The second cluster of solutions is the one of 
the "factorized" solutions, namely BHLUMI and SABSPV in its factorized 
versions (see eqs. @ and (||3))-EI All these solutions have the common feature 
that they include the bulk of the 0(a £ C) contributions in some form. It is 
worth noting that all of them lye in the band a re f ± 0.1%. The comparison 
between the results of eqs. (|) and ( [XQ|) is a measure of the degree of approxi- 
mation of eq. (|4]) with respect to the new improved solution of eq. ([10|) . The 



relative difference between them is maximum at z m i n = 0.1, where it is roughly 
0.1% and can be attributed to convolution effects, and minimum at z min = 0.9, 
where it amounts to around 0.02% since the convolution effects switch off nat- 
urally. These differences should be compared with the corresponding relative 
difference with the cluster of additive solutions, which on the contrary is larger 
than 0.1%. This means that the approximate solution of eq. (f|) is, after all, 



1 Actually, the results shown in Fig. [j] have been obtained by means of the sum of eq. (^) 
and eq. (||), which anyway differs from eq. ( flO| ) starting from 0(a 3 C e ) corrections; these 
corrections are numerically irrelevant, as checked for the case under consideration, and shown 
later on for a realistic ES. 



10 



a good approximation to eq. (|i~0D, which anyway supersedes it. It is worth 



stressing that the main difference between eq. (fD and (10) does not lye in 
the dynamical content, but rather in convolution and/or phase-space effects, 
that are correctly treated in ( |T0"|) while approximated in (f|). Note that the 
new improved solution and the reference cross section differ from one another 
by less than 0.03% over all the z min range. As far as the physical content is 
concerned, the results by NLLBHA are comparable with the factorized results 
by BHLUMI and SABSPV. In particular, also these results differ from the 
reference cross section by less than 0.1% over all the z m i n range and less than 
0.06% for realistic z min , namely 0.3 < z m i n < 0.7. The maximum difference be- 
tween NLLBHA and BHLUMI/SABSPV is at z min = 0.9, and could be traced 
back to the infrared- sensitive terms of eq. (0). This item deserves further 
investigation. 

The formulation described above in eq. (|T0"D is conceived for a Bare ES. 
Going to a Calorimetric ES requires a completely different treatment of final- 
state radiation in the LL part of the result. A first step can be simply drop- 
ping final-state radiation, which means putting f3f — in all the structure 
function corrections, consistently with the fact that, for an inclusive measure- 
ment, final-state mass "singularities" are absent. An improvement with respect 
to this choice can be taking into account final-state radiation having defined 
a proper final-state radiation factor. This can be done by considering that 
adding the correction due to photons collinear to final-state leptons results 
into a modification of the f3 factor [14] according to 

ft = 2(c0r) [ln(-t/m 2 ) - ll -+ ft = 2(c0r) In (^-) , (14) 



\5 2 s 

where 5 can be reasonably taken as the minimum aperture of the final-state 
cluster. The results obtained for a CAL02 ES (see @ for details) are shown 
in Fig. ^, where the same conventions as in Fig. [I] have been adopted. The 
only difference is that, in this case, the results by NLLBHA are absent and the 
numerical results relative to the present formulation have been produced first 
switching off final-state radiation (Fig. 2a), and then including its effect using 
in the structure functions ft = ft with ft given by ([14]) (Fig. 2b). As in the 
case of the BARE ES, two clusters of solutions can be recognized. The first 
group of the additive solutions, due to the lack of 0(a e C) contributions, differ 
from the reference cross section of about 0.10-0.15% for 0.1 < z m i n < 0.7 and 
of about 0.3-0.5% at the extreme value z min = 0.9, where the differences can 
be expected to be enhanced by infrared-sensitive terms. The second cluster of 
the "factorized" solutions (BHLUMI and SABSPV in its factorized versions), 
which include the bulk of the (D(a & £) corrections, is contained within a band 
<j re /±0.1%, both without and with final-state radiation. In particular, for the 
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case of no final-state radiation (Fig. 2a), we show also the numerical results 
corresponding to the sum of eq. @ with eq. (||), which, although in additive 
form, does contain 0(a & £) contributions. As can be seen, the difference be- 
tween the new factorized version of SABSPV of eq. ( |10"D and the result given 
by the sum eq. (||) with eq. is not appreciable, even at this level of pre- 
cision, since the two prescriptions differ for sub-leading terms starting from 
0{a 3 C £ ). The largest difference between the improved versions of SABSPV 
and the reference cross section is present at z m i n = 0.9 in Fig. 2a, and it has 
to be ascribed to neglecting final-state radiation corrections in SABSPV, for 
this case. This interpretation is confirmed by the situation shown in Fig. 2b, 
where the predictions by SABSPV include the effect of final-state radiation. 
Actually, the new improved solution and the reference cross section differ, in 
this case, by less than 0.02% over the full z min range, similarly to the case of 
the BARE ES. As a general comment, it is worth noting that, wherever dif- 
ferences between solutions are present at the extreme value z min = 0.9, these 
are due, as we explicitly shown independently of the ES, to differently treated 
infrared- sensit i ve cont r ibut ions . 

In conclusion, we have shown a new formulation for the calculation of 
photonic corrections to a generic kernel cross section. A method for computing 
the bulk of photonic 0(a e C) corrections has been proposed, and analytically 
tested versus theoretical results present in the literature for simple ES's. The 
algorithm has been applied to the precision calculation of the SABH cross 
section, relevant for the luminosity measurement at e + e~ colliders, and several 
numerical results have been shown and critically commented. 

The present formulation, as applied to the SABH process, supersedes the 
one developed in |7|, § and used in [Q, |3| as one of the tools for the estimate 
of the theoretical error on the SABH process itself. 

Following the strategy adopted in [§, |3|, an appreciable reduction of such 
an error will only be reached by means of a critical comparison of all the 
available formulations of the problem, for a wide ensemble of ES's. The present 
formulation is able to provide one of the ingredients for such an achievement, 
towards a theoretical error on luminosity well below 0.1%. 
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Figure 1: Comparison of MC and non-MC results. The relative differences 
between the codes involved in the comparison and the cross section by BH- 
LUMI taken as a reference cross section are shown as functions of the cut 
z min = E + E_/ El eam . E + _ are the energies deposited by the bare final state 
positron and electron, respectively. The details of the clustering algorithm 
(BARE1) are given in ||. The centre of mass energy is y/s = 92.3 GeV. 
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Figure 2: Comparison of Monte Carlo's. The relative differences between the 
codes involved in the comparison and the cross section by BHLUMI taken as a 
reference cross section are shown as functions of the cut z min = E + E_/ E% eam . 
E+- are the energies deposited in the positron and electron clusters, respec- 
tively. The details of the clustering algorithm (CAL02) are given in ||. The 
centre of mass energy is y/s = 92.3 GeV. 
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